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We analyze quantitatively the accuracy of eigenfunction and eigenvalue calculations in the frame 
work of WKB and instanton semiclassical methods. We show that to estimate the accuracy it 
(J_) ' is enough to compare two linearly independent (with the same quantum number) solutions to the 

^ ' Schrodinger equations with the potential under study, and with the approximating piecewise smooth 

potential. The main advantage of the approach is related to the appropriate choice of the approx- 
imating potential, providing absolutely convergent majorant series for the solutions. We test our 
method for a strongly anharmonic one dimensional potential, but the basic ideas inspiring our work 
and its results can be applied to a large variety of interesting chemical and physical problems which 
are of relevance to various molecular systems. 
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PACS numbers: 05.45.-a, 72.10.-d 



I. INTRODUCTION 



^ , It is a textbook wisdom that if the de Broglie wavelengths A of particles are small in comparison with the charac- 
' teristic space scales of a given problem, then the problem can be treated semiclassically. The commonly used WKB 
' method (phase integral approach) - Q is intended for the conditions of " geometrical optics" , in which the gradient 
of the action a is large, but slowly variable (this is suggested also by it containing the factor h, since we are dealing 
with a semiclassical approximation, in which h is taken as small). The corresponding condition can be formulated 
more quantitatively [3 as follows: A must vary only slightly over distances of the order of itself 



d{X(x)/2TT) 



dx 



« 1 , (1) 



I . where A(a;) = 2Trh/p{x), and p{x) is a classical momentum. However, this simple criterion is not a practical tool 
''O ' to estimate how accurate could be found semiclassical solutions of particular problems, since nothing is specified 
S . regarding the convergence of the semiclassical series. The criterion |^ does not work to estimate the magnitude of 
^ ' the error involved in the approximate calculation of physical quantities (e.g., matrix elements), neither to find the 
. . , domain of validity in the complex plane in which the semiclassical solutions are defined. Indeed from |^ one can 
^ ' conclude only that higher order corrections to semiclassical wave functions are small in the asymptotic regions, but 
this mathematical criterion has almost nothing to do with say the physical accuracy of semiclassical matrix elements 
which depends on the wave function accuracy in space regions providing main contributions into the matrix elements 
under consideration. For example, the energy eigenvalues are determined by the asymptotical region of the linear 
turning points (i.e., the region distant from these points), and as well by the proximity region to the second order 
turning points, since in the both regions the wave functions possess the largest values. 

Within the WKB method such kind of a physical accuracy estimation has been performed long ago by N. and 
P.O. Froman Q. They analyzed higher order corrections to the semiclassical wave functions and found that although 
those are really small over 1/7^ (7 ^ 1 is semiclassical parameter), the corrections are proportional to the factor 
\{E/^) — t7]~^, (where E is energy and U is potential), and thus the function has non-integrable singularity at the 
linear turning points where \{E/^) — t/] = (or, within the alternative to WKB semiclassical formalism so-called 
extreme tunnehng trajectory or instanton instanton approach - lUj , the corrections are singular in the second order 
turning points) . To surmount this problem in (see also , and |l4j ) the analytical continuation of the correction 
function into the complex plane has been proposed, and it gives impractically bulky expressions even for simple model 
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potentials. Since this problem has relevance far beyond WKB treatment of a particular model potential, this is an 
issue of general interest to develop a simple and convenient in practice quantitative method to study the accuracy of 
the semiclassical approach, and it is the immediate motivation of the present paper to develop a systematic procedure 
how to do it. 

The idea of our approach is to construct two linearly independent continuous (with continuous first derivatives) 
approximate solutions to the Schrodinger equation, which in the asymptotic region coincide with semiclassic solutions, 
and in the vicinity of the turning points - with the exact solutions of the so-called comparison equation (i.e. the exact 
solution of the Schrodinger equation for the chosen appropriately approximate near the turning points potentials 
Vc{X), henceforth will be referred to as the comparison potential). Although, scanning the literature we found one 
rather old paper [T^ with a similar comparison equation approach, but our accuracy criterion is formulated as the 
majorant inequalities for a certain matrix (which we find in the explicit analytical form and calculated numerically) 
connecting our approximate and exact solutions in the finite space interval (not only in the vicinity of the turning 
points). Since this method has largely gone unnoticed in the study of semiclassics, we found it worthwhile to present 
its derivation in a short and explicit form, and also to point out its practical usability. 

The remainder of this paper is organized as follows. In section ^1 we present the basic expressions necessary for 
our investigation. In this section we also present the main steps and qualitative idea of our method. Section IIIII 
contains our results. We derive the inequalities which enable us to find the finite space interval (not at the isolated 
points) where the solutions have to be matched, and calculate the 2x2 coordinate dependent matrix connecting the 
approximate and exact solutions. Since the semiclassical solutions of the harmonic potential coincide with the exact 
solutions, the accuracy of any semiclassically treated problem depends crucially on its potential energy anharmonicity. 
That is why as the touchstone to test our method the results presented in the section UTTl are applied to an anharmonic 
oscillator in section Hvl We end with some brief conclusions in the same section. 



II. SEMICLASSICAL EQUATIONS IN THE WKB AND INSTANTON FORMS 

Technically the basic idea how to overcome the difficulty of the semiclassical solutions in the vicinity of the turning 
points is reduced to an appropriate (admitting exact analytic solutions) approximation of the potential near the 
turning points. After that step one has to match the asymptotics of this exact solution to the Schrodinger equation 
for an approximate potential with the semiclassical solutions to the Schrodinger equation for the potential under 
consideration (i.e. approximate solutions of the exact potential) far from the turning points. To illustrate main 
ideas of any semiclassic method (and to retain compactness and transparency of expressions) we discuss here a one 
dimensional case. As it is well known jlj in the WKB method solutions to the Schrodinger equation are sought in the 
form 

^ = Aexp(^^), (2) 

where for the function a called action the one particle Schrodinger equation (traditionally termed as Hamilton Jacoby 
equation) reads as 



1 rda^ ^ 



2m V dx 



E-U, (3) 



where m is a particle mass, E is its energy, and U is external field potential. Since the system is supposed quasi- 
classical in its properties, we seek a in the form of a series expanded in powers of h. Depending on normalization 
prefactor A{x) entering Q can be also found but the corresponding equation (referred traditionally as transport 
equation) plays a pure passive role since it is fully determined by the action a found as the solution of the Hamilton 
- Jacoby equation 



ih 
m 



ISV , OA da 
_l_ 

2 dx'^ dx dx 



2m dx^ ' 



where in the spirit of the semiclassical approximation the last term (oc ti^) is neglected. Technically of course more 
convenient to use instead of h an expansion over equivalent but dimensionless parameter 7"^ <C 1 we will call in what 
follows as semiclassical parameter and define as 



(5) 
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where oq is a characteristic length of the problem, e.g. the tunneling distance, JIq is a characteristic frequency, e.g. 
the oscillation frequency around the potential minimum. Evidently the semiclassical parameter 7^1, and by its 
physical meaning it is determined by the ratio of the characteristic potential scale over the zero oscillation energy. 
We put h — I, and use f2o and oq to set corresponding dimensionless scales, i.e. we introduce dimensionless energy 
e = E/^D,o, dimensionless coordinate X = x/ao, dimensionless potential V = U/Qq (except where explicitely stated 
to the contrary and dimensions are necessary for understanding or numerical estimations). 

The analogous to |(2Jl, Q procedure for the Schrodinger equation in the imaginary time (instanton formalism, 
corresponding to the Wick rotation in the phase space, when coordinates remain real valued x x but conjugated 
momenta become imaginary — s- ipx) can be formulated as the following substitution for the wave function (cf. to 
©) 

tp = Ae{X) exp(-7cr£;) , (6) 

where the action as and we use the subscript E to denote so-called Euclidean action obtained from the WKB action 
a after the Wick rotation. Performed above rotation is not a harmless change of variables. The deep meaning of this 
transformation within the instanton approach is related to redistribution of different terms between the Hamilton - 
Jacoby and the transport equations. Indeed, like that is in the WKB method, eigenvalues for the ground and for the 
low-lying states are of the order of 7", while all other terms in Q are of the order of 7^. Therefore to perform a 
regular expansion over 7^^ for the substitution (jSJ one has to remove the energy term from the Hamilton - Jacoby 
equation, and to include this term into the transport equation. Besides in the first order over 7^^ one can neglect the 
term with the second derivative of the prefactor. As a result of this redistributions the both equations are presented 



2\dX J 

instead of the WKB Hamilton - Jacoby equation J^l, and the transport equation is 

OAe dcFE 1 d'^cjE 



V{X), (7) 



dX dX 2 dX^ 



Ae = eAE . (8) 



One can easily note by a simple inspection of the WKB Q , (QJ and of the instanton , (|Hl equations that although 
the both semiclassical methods can be formulated neglecting terms of the order of 7~^, therefore possessing the same 
accuracy over 7"^, the solutions evidently coincide in the asymptotic classically forbidden region V{X) ^ e/7, but 
their behavior, number and type of turning points (where any semiclassic approximation does not work) are quite 
different. For example in the WKB formalism there are two turning points where V{X) — (e/7) = around each 
minimum of the potential, while in the instanton approach since the energy does not enter the Hamilton - Jacoby 
equation (thus one can say that no classically accessible regions at all) the turning points are extremal points (minima 
for the case) of the potential. Furthermore as a consequence of this difference, in the WKB method all turning points 
are linear, whereas in the instanton approach they are second order (quadratic over X). 



III. ACCURACY OF SEMICLASSICAL APPROXIMATION 



Armed with this knowledge we are in the position now to construct our approximants. Let us introduce besides the 
comparison potential T4(X), one more specially chosen potential Vsc (henceforth will be referred to as the semiclassical 
potential). This potential is chosen by the requirement that the exact solutions to the Schrodinger equation with 
Vsc coincide asymptotically with the semiclassical solutions to the Schrodinger equation with the potential V{X) the 
problem under study. Thus according to the construction, the semiclassical wave function satisfies the equation 

and from here we can relate the semiclassical potential {Vsc{X)) with the bare one {y{X)) 

e«^nx)^^^-'(|^). (10) 

in the vicinity of the first order or of the second order (superscripts 1 or 2) turning points. 
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bmce near the turning point Xq the prefactors A^^) cx |X - Xq]-'^/^, and A^^) cx (X - Xq)" (where n is an integer 
number which occurs from the transport equation ((SJ solution at the energy e = n + (1/2)) the potential Vs^P at 
X — > Xq is singular and negative, and V"ic ^ has the same singularity (oc {X — Xq)^"^) but positive. The difference is 
due to the fact that near the WKB linear turning points we have deal with the vic ^ well, whereas near the second 

(2) 

order instanton turning points one has to treat the potential barrier Vsc ■ It might be useful to illustrate the essential 
features of the introduced above potentials 14 and Vsc applying the definition H10() to a simple (but the generic 
touchstone) example of the following anharmonic oscillator 

1/(X) = + + (11) 

We show in Fig. 1 the semiclassical and the comparison potentials associated with (IllH for the WKB (Fig. la) and 
instanton (Fig. lb) methods (a — —1.25, /3 — 0.5, and the energy window corresponds to n = 3 excited state of the 
potential dJ). 

The key elements to construct our approximants are the following combinations related to probability flows to and 
from the turning points 

^(^i) = ^'s} - = 27'*s-cH^i)*c-H^i) P ^s.{X)^c{X){ysc{X) - Vc{X))dX, (12) 

where Xi < Xq and analogously for X2 > Xq the flow function J{X2) is given by p2|l where the integration limits 
are from X2 to +00. Since the exact wave functions are continuous with continuous first derivatives (providing due 
to these features the continuity of the density probability currents), the idea of our procedure is to require the same 
from the approximate wave functions. 

The integrals entering J{Xi) and J(^2) can be calculated easily for any form of the potential, and the maximum 
accuracy of the any semiclassical approach can be achieved upon the matching of the approximate solutions at the 
characteristic points Xf2 where J(X*2) — 0. The points Xf'2 do exist in the case when the potentials Vc and Vsc 
intersect in the region where the approximate wave functions ^'^c and ^c are monotone ones. It is easy to realize 
(see e.g., Fig. 1) that the both points occur in the vicinity of the linear turning points for the potentials with 
d?V/dX^ > 0. One such a point disappears when the potential turning point becomes the inflection point, and there 
are no points at all for (PV/dX^ < 0. In the vicinity of the second order turning point the comparison potential 
Vc is a parabolic one. The curvature of the latter potential can be always chosen to guarantee the two intersection 
points always exist. The choice of the comparison potential corresponds to a certain renormalization (cx 7^^) of 
the characteristic oscillation frequency {d'^V/dX'^)x=Xo- Note in passing that the approach we are advocating here 
conceptually close (although not identical) to the scale transformation proposed by Miller and Good |0| and further 
developed infTij. 

Thus the conditions J{xf2) ~ allow us to construct well controlled approximate solutions to the Schrodinger 
equation. The accuracy of the approximation depends on the deviation of the approximate wave functions from the 
exact ones in the vicinity of the characteristic points Xf2- In own turn, the deviation is determined by the higher 
over {X — Xq) terms of the potential V{X) which are not included in the harmonic comparison potential Vc- Include 
explicitely the corresponding higher order terms to distinguish Vc and Vsc potentials, we find in the vicinity of the 
linear turning points 

yii) -Vc^-^[X- Xo)-' + ^ (X - Xof , (13) 
7^ 2 

where the universal numerical constant ci = 5/32, and the second term in the r.h.s. is related to deviation of the bare 
potential from the linear one. The same manner near the second order turning points 

VjP -Vc~-^{X- XoT^ + a{X - Xof , (14) 
T 

where the universal constant C2 — n{n — l)/4 is zero for the lowest vibrational states n — 0, 1, and the last term 
describes non-parabolicity of the potential. Note that unlike the semiclassical action which within the instanton 
method is independent of quantum numbers n, the position of the characteristic points does depend on n, and the 
Xf2 points are placed near the boundaries of the classically accessible region. 
Now we are in the position to construct the approximant wave functions 
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which are the solutions to the Schrodinger equation with the following piecewise smooth approximating potential 

^^''^-\Vs..x<x*^x>x* ■ ^''^ 

The wave functions calculated according to (|15|1 in the framework of the instanton approach close to the Weber 
functions in the classically accessible regions, but their exponentially decaying tails in the classically forbidden regions 
correspond to the exact (bare) potential, not to its harmonic approximant. Analogously in the WKB method these 
functions IjlSI) coincide with the semiclassical ones out of the interval {xf , xf)^ and with the Airy functions in this 
interval. 

To proceed further on we have to relate our approximant wave functions IjlSfl and two linearly independent solutions 
to the bare Schrodinger equation 'Si and ^'2. It can be written down formally as 

^{X) = ^{X)+ [ dX'v{X')G{X,X')'^{X'), (17) 

where v — V{X) — V{X), G{X,X') is the Green function for the Schrodinger equation with the potential H16|l. 

GiX,Xi) = const[^i{X)^2{Xi) ~ ^i{Xi)^2{X)]; Xi < X , (18) 

where the constant in H18() is the Wronskian equal to (27)^^ in the instanton method, and 1(27)"^ within the WKB 
approach. In (|17|l 5* = (^'i , ^1^2) (the same definition for Xf), and we take the turning point Xq, where the functions 
^ and ^ are close to each other at the lower integration limit. 

The solution to the integral equation (|17|l is expressed as the Neumann series expansion, 



■^{X) = ^{X)+ [ dXiv{Xi)G{X, Xi)*(Ai) + • • • , (19) 

JXn 



<!-[ I dX^v{X,)G{X,X^)^{X^)\ . (20) 



rX 
IX, 

and the m-th order term can be factorized and estimated as 

r-X 
IXo 

The integrals entering this estimation 



TO 



L 



X 

dXiv{Xi)G{X,Xi)^iXi) = Li2/22(X)*i(^) + LU/21{X)^2{X) (21) 



IXo 

contain the 2x2 matrix with the following matrix elements 

U,= [ dX'^,iX')v{X')^,iX'). (22) 

JXo 

It is convenient to introduce the matrix C*^"^ relating the n-th order wave function correction with the wave 

function \I' 

'««Ucwf|.V (23) 



and the full connection matrix between the exact and approximate wave functions 
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According to the inequality (|20(l 

Cl") < 1q , (26) 



where 

Co = ( ) . (27) 

\ i^22 ^12 J 

Combining finally the expressions (|25() - H27|l we end up with the upper and lower bounds for the correction matrix 
C estimation 

l + Co<C < exp(C'o) . (28) 

We conclude from (|20|1 . H22|l that the Neumann series posses an absolute convergence if all the matrix elements Lij 
are finite. Besides, unlike the correction functions introduced within the Froman approach the integrals in H22(l 
have no singularities on the real axis. Evidently the integrals H22II are finite with the oscillating WKB functions, since 
the perturbation potential v{X) is not zero only in the close proximity to the characteristic points Xf'2- However in 
the instanton method due to mixing of increasing and decreasing exponents, the matrix elements L22 is divergent. 
Despite of this divergency the product L22VE'i, we are only interested in, is finite, and it is convenient to perform one 
more transformation to exclude explicitely this divergency. 

Technically one can easily eliminate the both off-diagonal elements of the matrix Cq and thus to get rid of the 
divergency of the (exponentially decreasing solution ^i) amplitude due to the contribution to the the exponentially 
increasing solution ^'2. These linear transformations renormalize the correction matrix elements L22 and Ln as follows 

L22{X)^i{X) = L;2{X)-^2{X) ■ L^i{X)^2{X) = L\^{X)^i{X) , (29) 

where the renormalized matrix elements L22 and L\i read as 

L*22{X) = dX'^i^^^^^^*i(X')w(X')^2(X') < Li2{X) , (30) 

' Jxo *i(X')*2(^) IV ; V ; 2V 7 - V 7 

and 

Ll.iX) = dX'l^^^^p^^^iX')viX')^2{X') > L,2iX) . (31) 

JXo Wi(A)W2(A') 

Now all the integrals entering L|]^, L22, and L12 are convergent and finite at any X, the correction matrix Cq is 
transformed into the diagonal and positively defined matrix Cq 

The n-th order corrections in the instanton approach satisfy the inequality 

< S¥"^ < {niy^Ln - ii2)*i , < (J^a"^ < {niy^L^ - i22)^'2 • (33) 

Explicit summation of r.h.s in (|33|l gives us the upper and the lower bound limits for the solutions of the initial 
Schrodinger equation, i.e. the stripe where increasing and decreasing solutions are confined 

|*i(^)| < |*i(^)| < |*i(X)| exp(LJi - L12) , \^2{X)\ < \^2{X)\ < 1*2(^)1 exp(Li2 - L*^) . (34) 

It is our main result in this paper, and the stripe 134|) gives the accuracy of the semiclassical instanton method. Besides 
we are in the position now to estimate the contribution of increasing semiclassical solutions into decreasing ones (what 
is relevant to solve eigenvalue problems). The summation convergent majorant series enables us to estimate the upper 
bound for this contribution 

§i(l-exp(-ii2)). (35) 
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Therefore at L12 <C 1 the summation of all order perturbation terms enhances the 1-st order correction by the factor 
^12/^22- Analogously the majorant estimates described above can be used to construct the connection matrices 
linkinff the semiclassical solutions through the turning points. The comparison of the bare connection matrices (see 
e.g., Q, 0, 01) with the matrices calculated accordingly to H22I) - l|3U|l provides the estimates for the eigenvalue 
accuracy. 

As it was mentioned already, for the WKB method the procedure is even more simple, since no any divergency and 
therefore no need to perform the transformation (|29|l . The similar to (|20|) - (|21|l factorization gives the Cq matrix (cf. 
with H27|l for the instanton approach) 



Co = i 



L12 —L 



11 



L22 —L12 

and the estimations for the n-th order contribution (cf. with H33|) can be formulated now as 

1 /AO" 



1(7(2") I < 



(2n)! V A 



and 



|(7(--)|<^^A"Co 



(36) 



(37) 



(38) 



where we denote A = ~ ^11-^22- Correspondingly to (|37|l . (|38|l the diagonal and off-diagonal correction matrix 
elements are bounded from above 



|Cll/: 



22 



< 



/ t 

cos V A + iLi2- 



and 



1^*12/211 < 



ill 



sin \/A 



A 



(39) 



(40) 



The whole procedure we employed is rationalized in the Fig. 2, where we compare the solutions to the compari- 
son equation with the anharmonic oscillator semiclassical wave functions computed within the instanton and WKB 
approaches and indicate the optimal matching points X'^ found accordingly to the condition J{X'^) = 0. 



IV. ANHARMONIC OSCILLATOR 



In closing let us illustrate how our estimations (|34|l . (|38|l work for a strongly anharmonic potential Although 
it is not great triumph to re-derive the known results, our derivation illustrates several characteristic features of 
the correction matrix techniques derived in the section IIIII better accuracy, rapid convergence, simple disposal of 
divergences, and ease of computation in particular. The main message of our consideration in the precedent section 
IIIII is that the quantitative accuracy of the scmiclassics depends crucially on the proximity of the semiclassical wave 
functions to the solutions of the comparison equation in the region of the asymptotically smooth matching. Therefore, 
it is tempting to improve the accuracy by taking into account the anharmonic corrections to the comparison potential 
Vc{X). However, since the eigenvalues and the normalization of the wave functions are almost independent of the 
detailed behavior in the vicinity of the linear turning points (because near these points, situated at the boundaries 
of the classically accessible region, the probability density (i.e. jvE'P) is exponentially small) this idea is useless for 
the WKB approach. In contrast with this, for the instanton method, the accuracy can be improved considerably 
upon including the anharmonic corrections into the comparison potential. Indeed, within the instanton method the 
accuracy is determined by the vicinity of the second order turning points where the wave functions acquire the largest 
values (and just in this region the smooth matching described above has to be performed). 

Let us remind first the traditional (but formulated within the semiclassical framework) perturbation theory. Keeping 
in the transport equation the second derivative of the prefactor A, the system of the equations of Hamilton - Jacoby 
(0) and the exact transport equation 



dX^ 



27 



dA da 
dXdX 



d^a 
dX^ 



-e\A 



= 



(41) 
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are exactly equivalent to the Schrodinger equation under consideration. For the second order turning points the 
anharmonic corrections 

Vp = \ {aX^ + PX^) (42) 

can be considered as a perturbation and it is convenient to include this perturbation Vp into the transport equation 
(I41|l . Then the comparison equation is reduced to the inhomogeneous Weber equation, and its solutions can be 
expanded over the Weber functions D^{X) [23| (see also ^13], 14]) 

*,(X) = [du{X) + ^ hkuD,+k{X)^ , (43) 

where N~'^ = 1 + is the wave function normalization factor, the expansion coefficients are proportional to the 

small parameters a/^, and P/j, and the Weber function index f is related to the energy eigenvalue e = + (1/2). 
This expansion (|43|) looks like a conventional perturbation series, but it does not. In the comparison equation we are 
keeping the both (decreasing and increasing) waves, and as a result of it, the indices of the Weber functions v -\- k are 
not integer numbers. In the first order over the perturbation Vp the only non-zero coefficients in (|43() correspond to 
the following selection rules 

fc = ±1,±3; andfc = 0,±2,±4 (44) 

for the cubic and fourth order anharmonic corrections respectively. Explicitely these non-zero expansion coefficients 
can be found by straitforward calculations, and they are 



6-3,. = -^au[v - - 2)(!. - 3) , 63. 6-4. = 2){y - - 4) , 64. = . 

On equal footing we can find the perturbative corrections to the Bohr-Sommerfeld quantization rules, and therefore 
the eigenvalues. The calculation is straightforward, though deserves some precaution and rather tedious. Skipping a 
large amount of tedious algebra we end up with the fractional part of the quantum number v 

t^^n + Xn, (46) 
and up to the second order over the anharmonic perturbation Vp we find 

However the described standard perturbative approach leads to qualitatively wrong features of the solutions. For 
example, the wave functions l|43|) , l|45|l are represented as a product of v independent exponential factors and dependent 
of u polynomials. As it is well known in one dimension the n-th excited state wave function must have n zeros 
(and the number of zeros may not be changed by any perturbation). However in the m-th order perturbation theory 
approximation, the wave function (|43|l corresponding to a certain excited state n contains Hermitian polynomials up 
to the order n -\- 3m or n -I- 4m for the cubic or quartic anharmonic perturbations, respectively. Therefore some false 
zeros of the wave function appears in the standard perturbation theory, and the region where the function oscillates 
becomes more and more wide in the higher order over perturbations approximation. The contributions of these 
qualitatively and quantitatively incorrect higher order terms become dominating in the asymptotic region at \a"f \ 1, 
/37 ~ 1. It conforms with the classical results due to Bender and Wu |0| who have shown that for the quartic 
anharmonic potential (a = 0, /3 > in (|ll|l ) the convergency radius is zero. Moreover, for /S — and for an arbitrary 
small a (|llll is the cubic anharmonic potential, i.e., the decay one. Thus it should have only complex eigenvalues, 
what is not the case for the eigenvalues calculated within the perturbation theory. The method we developed in 
section iHll enables us not only to estimate more accurate the anharmonic corrections to the eigenvalues, and to bring 
the whole schema of the calculations in a more elegant form. Our finding of the correction matrices is not merely to 
surpass a technical difficulty of the standard perturbative method, it is more one of principle, and we will show that 
the method has no drawbacks of the perturbation theory. 
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The proof proceeds as follows. Let us consider first the instanton method for the anharmonic potential H41() 
possessing one second order turning point X = 0. As it was shown in the section UTTI one has to find also two other 
characteristic points which are the roots of the equation (|12|) . We denote the points as xf, and (to refer by the 
self-explanatory subscripts L and R to the left and to the right from the turning point X = 0). At the next step using 
the correction matrices introduced in the section Hill we can define formally the transformation of our approximate 
wave functions 1)15(1 into the unknown exact wave functions 5* as 



where as above the subscripts L and R refer to the wave functions in the regions to the left and to the right from 
the turning point X — 0. We do not know the correction matrix C but we do know (see H23|l - ((31|l l the boundary 
estimations for the matrix. 

The Eq. (|48|) can be used also to correct the known at the second order turning point the connection matrix M 
Indeed the connection matrices link the semiclassical solutions in the AT-regions to the left and to the right 
from the turning points. For the isolated second order turning point (we are dealing within the instanton method), 
the connection matrix M links the exponentially increasing and decreasing solutions in the space regions separated 
by the turning point. The condition ensuring the correct asymptotic behavior is the quantization rule for this case 
which can be formulated as Mn = {Mij are the matrix elements of the connection matrix M). Since in the regions 
to the left and to the right from the turning point our approximate solutions ^' coincide by their definition ((15(1 with 
the semiclassical ones, the correction matrix method enables us to correct the quantization rule too. Namely, the 
quantization rules are formulated within the connection matrix technique read now as 

Ci^T2C^2 - C^iC^,^-^^^^^ + {Cgci + C^.Cg) cosin:.) = , (49) 
where the Stokes constant for the second order turning point Q is 



and C^'^ are the correction matrices at the AT]^ or xf^ characteristic points respectively. Expanding the Gamma 
function entering ((49(1 around the integer numbers, i.e., as above ((45(1 . p = n + Xn one can find from the equation 
((49(1 the fractional part of the quantum number. If we were known the correction matrix C the solution of l(49() would 
provide the exact eigenvalues. But we do know only the estimations from below and from above for the C matrix. In 
the same spirit we can calculate the estimations for the correction matrices (and therefore for the eigenvalues) within 
the WKB approach. 

Luckily it turns out that the mathematical nature of the semiclassical problem is on our side here, and, in fact, even 
the first order estimation from below C^*'-' II23() gives already the accuracy comparable with the standard perturbation 
procedure, and the estimation from above ((25(1 gives the eigenvalues almost indistinguishable from the "exact" ones 
obtained by the numerical diagonalization of the Hamiltonian. The same true for the wave functions found by the 
correction matrix technique. We show in Fig. 3 l^'sl" for the same anharmonic potential ((11(1 with a = —1.25, P = 0.5. 
Clearly the exact numerical results and those obtained by our correction matrix techniques are correct qualitatively 
and in the very good quantitative agreement (indistinguishable starting from the second order approximation) unlike 
the situation with the standard perturbation theory. Besides we present in the table the eigenvalues of the anharmonic 
potential. We take the anharmonic coefficients a and (3 in 1(1111 so large, that corresponding perturbations of the eigen 
values are of the order of the bare harmonic frequency (one in our dimensionless units a = —1.2, /3 = 0.5). In the 
table the eigenvalues found by the numerical diagonalization are presented in the column /. The column // contains 
the harmonic approximation results, the column /// is the second order perturbation theory ((47(1 . and the columns 
IV and V results are obtained by applying our correction matrix technique: estimation from below with the first 
order correction matrix ((23(1 in the column IV , and the estimation from above with the matrix ((25(1 in the column V. 

To conclude, as we have shown how to estimate the corrections to the main technical tool for the semiclassical 
approach, the connection matrices linking the solutions to the left and to the right from the turning points. Everything 
(e.g., the upper and the lower bounds for Xn) is determined by the matrices L12 and L22. Approximating to 1(11(1 
potential is found from ((13(1 and after that straightforward computing according to ((2211 the matrices Ly and their 
renormalization ((30(1 leads to the corrections we are looking for, presented in the Fig. 4, which allow us to estimate the 
accuracy of the semiclassical eigenstates and eigenf unctions. We conclude that even for a strongly anharmonic potential 
the both methods (WKB and instanton) are fairly accurate ones (about 5%) up to the energy close to the potential 
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barrier top (in the region of negative curvature, we already discussed above). It is worth noting that in the frame work 
of the conventional perturbation theory (due to zero convergency radius with respect to (3 coefficient in (|11|) ') pure 
computational problems to get the same accuracy become nearly unsurmountable, see e.g., 0, Our findings 

show that to estimate quantitatively the semiclassical accuracy it is enough to compare two linearly independent (with 
the same quantum number) solutions of the initial potential under study, and of the approximating piecewise smooth 
potential. The main advantage of the approach is related to the appropriate H13(l , (|16|l choice of the approximating 
potential, providing absolutely convergent majorant series (|17|) for the solutions. Actually our correction matrix 
technique is a fairly universal one and enables to estimate (and improve!) the semiclassical accuracy for arbitrary one 
dimensional potentials with any combination of the turning and of the crossing points. 
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Figure Caption 

Fig. 1 

The characteristic semiclassical potentials T4c (dot-dashed hnes) and Vc (dashed hnes) for the bare anharmonic 
potential a = -1.25, f3 = 0.5 (7 = 33, and the energy window corresponds to n — 3): 

(a) instanton approach; 

(b) WKB method. 
Fig. 2 

Semiclassical wave function 'i's {n — 3) for the anharmonic potential (|ll|l with a = —1.25, j3 — 0.5 (7 = 33). Stars 
indicate the matching points, dashed lines show the solutions to the comparison equations, and: 

(a) - solid line traces the instanton solution; 

(b) - solid line shows the WKB wave function. 
Fig. 3 

Comparison of the exact l^'ap (solid line) for the anharmonic potential Hll|l with a — —1.25, (3 — 0.5 (7 = 33) with 
two lowest (zero and first order) approximations of the correction matrix method (dashed and dot-dashed lines). Note 
that the second order approximation with the relative accuracy 10^^ is indistinguishable from the exact numerical 
results. 

Fig. 4 

Corrections to the decreasing solutions in for the anharmonic potential (|ll|l : 
/3 = 0.5 , n = 3 , 7 = 33; 

(1, 3) - instanton method, (2, 4) - WKB, (1, 2) - the first order corrections Ln, (3,4)- the upper bound estimation 
summing up all terms. The vertical dashed lines indicate the values of the cubic anharmonic term (a) where the 
inflection point and new extrema of the potential are appeared. 

Table 

Eigenvalues of the anharmonic potential (|11|) {a = —1.2 , (3 = 0.5). 
/ the eigenvalues found by the numerical diagonalization; 
// the harmonic oscillator eigenvalues; 

/// the eigenvalues in the second order perturbation theory H47I) : 

IV the eigenvalues estimated from below by the correction matrix H23|l: 

V the estimation from above with the matrix H25|) . 
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